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Abstract 

We investigate multidimensional model for incompressible 
non-Newtonian fluids. Using method of energy estimates we prove 
the property of finite speed of propagations of the solution support for 
this problem. We find sharp bounds of the propagations by L 2 -norm 
and L 1 -norm of initial data. 
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1 Formulation of the problem and main results 

There are many substances like geological materials, liquid foams, polymeric 
fluids, which are capable of flowing but which exhibit flow characteristics 
that cannot be adequately described by the classical linearly viscous fluid 
model. In order to describe some of the departures from Newtonian behav- 
ior evinced by such materials, many idealized material models have been 
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suggested. There are several phenomena which appear when studying non- 
Newtonian fluids: the shear thinning and the shear thickening, the ability 
of a creep, the ability to relax stresses, the presence of normal stress differ- 
ences in simple shear flow, the presence of yield stress. For more details see 
[T2] . The equations governing the unsteady motion these fluids have been 
previously investigated, e.g. (5j QUI CD] , etc. 

Let components of the symmetric deformation velocity tensor are given 

by 

1 f dui i duj 



(| Du | = (Duij&Uij)^). 

We consider multi-dimensional problem of non-Newtonian model which 
has the following form 

{u t + (u ■ V) u = ^ l div{\Du\ p - 2 T)u) + Vtt, (1.2) 
divu = 0, (1.3) 
u(0,x) = u (x). (1.4) 
Here u is the velocity field, n is the pressure for an incompressible power-law 
fluids, u is the initial value of the velocity. Let Hi > 0, 

M e (W 1 * 2 ^))* n H, supp u , k (x) Q R N , k = L~iV, (1.5) 

where := {x = (x',x N ) G R N : x N < 0} and H = {u G L 2 (R N ) N : 
divu = 0}. 

Existence of a weak solution for p > and the uniqueness and reg- 
ularity for p ^ ^g^, iV = 2,3 were proved in [5], [8]. The more detailed 
information about qualitative behavior of the solution can be found in [T]. 

Definition 1.1. The function u(x, t) = (ui(t, x), . . . , u^it, x)) such that 

ii G L p (0, T; W 1,P (R N ) N ) n L 2 (0, T; W^ 2 (R N ) N ) fl C(0, T; tf), 
u t eL 2 {0,T;H) 

is called a weak solution to problem (C) if for a.e. t > 0, the integral identity 
J u t (t)<fdx+ j \u-V)u<pdx + Hx J \Bu(t)\ p - 2 ~Du(t) :B^dx = (1.6) 

R N R N R N 

is satisfied for every if G L p (0, T; W 1,P (R N ) N ) n L 2 (0, T; W 1 > 2 (R N ) N ) with 
divip = 0. 
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Basic definitions and properties of p-Laplacian we can find in the books 
of Vazquez [T9], EH]. The evolution p-Laplacian is one of the most widely 
researched equations in the class of nonlinear degenerate parabolic equations, 
already studied by Raviart [13]. Large time behaviour of the solutions was 
investigated in Kamin, Vazquez |6] and they proved that explicit solutions 
found by Barenblatt in 1952 are essentially the only positive solutions to the 
Cauchy problem with initial data 

u(x,0) = MS(x), M > 0. 

They established that any nonnegative solution with globally integrable ini- 
tial values is asymptotically equal to the Barenblatt solution as t —>■ oo. In 
work of Lee, Petrosyan and Vazquez [7] the property of asymptotic concavity 
was proved. 

Now we give brief explanations of our method of proof of the finite speed 
of propagations. This method we call as the method of nonhomogeneous 
functional inequalities and it is some adaptation the energy method. It is 
connected with nonhomogeneous variants of Stampacchia lemma, in fact, it 
is an adaptation to higher order equations of local energy or Saint-Venant 
principle like estimates method. It was elaborated with respect to higher 
order quasilinear parabolic equations of the monotone type in [HI CG2 fT6] . 
Also it was developed with regard to thin-film equations with nonlinear con- 
vection in [TTj, HE]. Fluid mechanics is one of the most natural fields for the 
application of energy methods (see, [Tj El [3]). This is because the fundamen- 
tal conservation and balance laws employed for its description suggest an 
adequate choice of energy functions which, in turn, produces the formation 
of a free boundary. "Energy methods are special interest in those situations 
in which traditional methods based a comparison principles have failed. This 
method yields the formation of a free boundary; in other words, this means 
that the support of solution is localized in space-time domain." (see, pQ) 

The main results are the following. 

Theorem 1.1. Let u(x,t) be a weak solution of problem (C). Let p ^ 3 j ^ h 2 2 . 
Then there exists the function T(t) G C[0,T], T(0) = such that 

2 2p+N{p-3) 

T(t) = a max{t 2 p+ JV (p- 2 ) , t2p+iv(p-2)} \/t > (1.7) 

and 

supp u k (t, .) n {x = (x, x N ) G R N : x N ^ T(t)} = 0, k = T^N, (1.8) 
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where C\ = Ci(p, /xi, ||iio||x, 2 (K iV ) JV ) ^ s a positive constant depending on known 
parameters only. 



Remark 1.1. If (Tjjty is valid then a solution to (C) has the property of 
finite speed of propagation of the solution support. 

Remark 1.2. The estimate jX 71) is sharp for scalar equation if the initial 
function Uq(x) belongs L2($i N ) N (see, for example, [Lflj ). 



Theorem 1.2. Letu(x,t) be a weak solution of problem (C) . Let \\u\\l 1 (^n^n ^ 
IM|li(R*-)* and p is ^y+T • ^ en ^ e function T(t) has the following form 

1 p+iV(p-3) 

r(t) = c 2 max{tf+^- 2 ), tp+^(p-2) } Vt > 0, (1.9) 
where c 2 = C2(p,fix, ||mo||li(]r jv ) jv ) ^ s positive constant. 



Remark 1.3. 77ie exponent 1/ (p+N(p— 2)) /rom ( tJ.^I) is well-known Baren- 
blatt exponent for p-Laplacian scalar equation (see, for example, L^j). 



2 Proof of Theorem 1.1 



Lemma 2.1. Let u(x,t) be an arbitrary weak solution of problem (C). Let 

SjV+2 
N+2 



p ^ 3 ^to ■ Then the following estimate is valid 



sup / \u(t)\ 2 dx + — II \u\ 2 dx dt + jii II \Vu\ p dx dt ^ 



te(o,T) 



T 

fl(s+S) Qt(s+8) Qt(s+S) 



\u\ p dxdt+^ 1 1 \u\ 3 dxdt, (2.1) 



Qt(s) Qt(s) 

Vs > s > 0, 5 > 0, T > 0. 
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Proof. For an arbitrary sel 1 and 5 > we consider the families of sets 

Q(s) = {x = (x',x N ) e R N : x N ^ s}, Q T (s) = (0,T) x Q(s), 
K(s, 5) = Q(s)\Q(s + 5), #x(s, 5) = (0, T) x if(s, 5). 

Next we introduce our main sequence cut-off functions r) Sj s(x) G C 1 (R JV ), 
which possess the following properties: 

^ 7^(2; < 1 Vx e ffi , ^(x = I \ ; (2.2) 

|Vt7 S)5 K^ VxeK( S) 5). (2.3) 
o 

Test integral identity (11.61) by 

ip(x, t) = u{t, x)ip{t } x), (2.4) 

where u(t, x) is a solution of the problem (C) and, without loss of generality, 
we suppose that (V x ip(t, x), u(t, x)) = 0. Here 

ip{t, x) = T]sA x ) ex P ■ T_1 ) VT > 0, 

Then we obtain 
\ J J^W^ + AH J |Dn(t)| P V dx+ 

in J \Du(t) \ p ~ 2 (Du(t) Vip) u{t) dx + j \u- V)u(t) u{t)^dx = 0. (2.5) 

We consider integral terms in the last identity. After using Young's inequality 
we get 



j \Bu\ p - 2 (DuVij)udx ^e^xp^t-T- 1 ) J \Du\ p dx+ 

K(s,S) 

exp(-t-T" 1 ) J \u\ p dx, (2.6) 



K(s,S) 

c(gQ 

5* 

K(s,S) 
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J{u- V)uu(t)ipdx = iexp(-t-T _1 ) / \u{t)\ 2 u{t) Vipdx < 

exp (-t ■ T~ l ) J \u\ 3 dx, (2.7) 

Integrating (12.5ft respect with time, and taking into account (12.61) and 02.71) . 
after simple computation, we deduce that 

sup / |w(t)| 2 <ir+^ // \u\ 2 dxdt + 2fi 1 [[ \Bu\ p dxdt^ 
te(o,T) </ T ii J J 

Q(s+<5) Qt(s+S) Qt(s+S) 



\u (x)\ 2 dx+^r^- J J \u\ p dxdt+^ J J \u\ 3 dx dt+fii6i J J \Du\ p dxdt, 

n(s) k t {s,s) k t (s,s) k t (s,s) 

(2.8) 

where sGM 1 , 5 > 0, T>0. It follows from CLE) that 

J \u (x)\ 2 dx = Vs ^ s ^ 0. (2.9) 

f2(s) 

Finally, from (12^81) . (T2T9D we obtain 

sup / \u(t)\ 2 dx+l- ii \u\ 2 dxdt + 2fi 1 [[ \Bu\ p dxdt^ 
te(o,T) J TJJ J J 

n(s+S) Qt{s+S) Qt(sU) 

<{ ' l) ^ \u\ p dxdt+- [[ \u\ 3 dxdt + fuE! [[ \Bu\ p dxdt, (2.10) 



$ p J J 6 _ 

K T (s,S) K t (s,S) K t (s,6) 

Vs ^ So ^ 0, 5 > 0, T > 0. Choosing e\ > sufficiently small, and iterating 

the limit inequality (12101 ). we get (CD]). □ 

Proof of Theorem li.il We denote 

-Rt(s, 5) := — |w| p <ic <it + - J J \u\ 3 dx dt. 

Qt(s) Qt(s) 
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We introduce the energy functions related to our solution: 



A r (s) :-- 



\u\ p dxdt, Bt(s) :- 



\u\ dxdt. 



Applying Nirenberg-Gagliardo's (see, Lemma A. 2 of Appendix A) and Holder's 
inequalities, we get 



A T (s + 5) < cT^(R T (s,8))^\ a, = (3, - 2p+N[p _ 2) , 



P(p-2) 



B T (s + 8) ^ cT a \R T { Sl 8)) l+p \ a 2 



2p+N(p-3) a = y 
2p+JV(p-2)> > J2 2p+N(p-2)' 



Next we define the functions 

C T (s) := (A T (s)) 1+ ^ + (B T (s)) 1+ ^ 

Then 



where 



C T (s + 8)^ cF{T) {8- pP C^{s) + 8- p C^{s)), 



P = (1 + + F(T) = max{T ai(1+/32) , T aa(1+/3l) }. 



(2.11) 



Now we choose the parameter 5 > which was arbitrary up to now. First, 
we introduce the notations 



8^'is) 



2cF(T)Cfr{s 



p " , 5 { t\s) := 2~cF(T)C^(s) 

J T {s) := max{5!f\s),5P{s)}. 

We obtain the following main functional relation for the functions Jt(s) 

J T {s + J T {s)) ^ e J T {s) Vs ^ s > 0, < e < 1. (2.12) 

We set s = —28, 8 = s' > in (12.81) and pass to the limit as s' — > +oo. Using 
the boundedness of functions At(s) and Bt(s), we get 



C T (0)^cF(T)\\u \\l. 



(2.13) 



Now we apply Lemma A.l to the function Jt(s) of (12 . 121) . As a result, we 
get 

j T ( s ) =OVOs + ^r(so). (2.14) 
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Let s = 0. Then, in view of ( 12.1 3D . we find 

1 1 2 2p+iV(p-3) 

J T (0) < c max-fi^T)^ 1 ^, F(T) 1+/3 i } ^ c max-fT 2 ^?- 2 ) , T^+^p- 2 ) } 
VT > 0, c = c(p, /ii, iV, IKIk)- Choosing in (12T4D 

2 2p+iV(p-3) 

s = r(T) = c max{T 2 p+^(p- 2 > , T^+^tf- 2 ) }, 

where p > 3gM we obtain that J T (r(T)) = 0. Thus u k {T,x) = for all 
x E {x = (x',xn) '■ xn ^ r(T)}, /c = l,iV. And Theorem 1.1 is proved 
completely. □ 



3 Proof of Theorem 1.2 

Lemma 3.1. Xet u(x, t) be an arbitrary weak solution of problem (C). Let 
IMIliOr^)^ ^ II^oIUiOr^)^ and p ^ 3 ^+i ■ Then the following estimates for 
the decay rate are valid 

A T {s) + B T {s) «C c(0) T (s-^p- 1 ) + s -p+^-3)) (3.1) 

Vs > 0, T > 0, where c(9) = c(p, N, 6), 6 = ||wo|| Li ( M jv)jv- 

Proof. Applying the interpolation inequality of Lemma A. 2 in the domain 
K(s,5) to the function v = \u\ for a = d = p, b = 1, dz = c8~ N ^~ x ^ p , and 
integrating the result with respect to time from to T, we obtain 

u\ p dxdt s$ c5' N{p - 1) Te p + cT 1 ~ ei Q p{1 - dl) x 



K T (s,S) 

X ( II \ Vu \ Vdxd *) ' where ^i = ^S)- (3-2) 

Vt(s,<5) ' 

Similarly, applying the interpolation inequality of Lemma A. 2 in the domain 

K(s,5) to the function v = \u\ for a = 3, d = p, b = 1, <i 2 = c5 _2Ar / 3 , and 
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integrating the result with respect to time, we find that 

302 



\u 



dxdt^c 5~ 2N T 3 + c T 1 v x 



K T (s,8) 

3 



x\JJ \Vu\*dxdt j , where 8 2 = 3{p+ ^_ 1)y V > ^r v (3.3) 

K T (s,8) 

Inserting (13.21) and (13.31) in f!2 . lOfl and applying Young's '^"-inequality, we 
get 

L T (s+6) := sup / u 2 (t)dx+l- \\ u 2 dxdt+^C // \Vu\ p dxdt ^ 
te(o,T) i J Ji 

n( s +s) Q T (s+s) Qt(s+s) 

£4^1 // I Vm| p rfx c?t + £3 sup / |u(t)| 2 fix + // \u\ 2 dxdt+ 



*e(o,T) 



Kt(s,<5) -^(s,<5) K t (s,S) 



C £4 (6)T(r (p+iV(p - 1)) + 5-S^T ), (3.4) 



Ve 4 > 0, s ^ 0, 5 > 0, where C £4 (0) is a constant depending on p, N, e 4 
and 0. Choosing £j > 0, i = 3,4 sufficiently small, by the standard iteration 
procedure we establish that 

L T (s + 8 ) ^ C(&)TU(5 ) V s > 0, 5 > 0, (3.5) 

where 

p+iV(p-l) 

U(S ) := 5- (p+iV(? '- 1)) + 5 Q i**fc=tf 

Let 5 -> +00 in (Q, and use ([3J)]) for s = and 5 Q = s > 0. We 

eventually obtain 

A T (s) ^ C{Q)TU dl {s), B T {s) < C(0)T?7 3e2 / p (s), 

for every s > 0. Hence (13. ip follows. 

□ 
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Proof of Theorem M.2. From (12.141) and the decay estimations (13.11) , we get 

G(s ) := so+j^Mso) ^ s +c m a x{F^(T)(A™(s )+B™(s )), 
F^(T)(^ +/51 (s ) + B± +02 (s ))} + c(6) max{F^(T)T^Ws p{1+M , 

i & mil i ^ ^fe(p-i) 

F?(T)T^+fes tp+iV(p - 3)K1+/J2) } < G(s ) := s + 

c(6) max{T^Ws p(1+/5l) , T ^Ws (p+JV(p ~ 3))(1+,32) }. 
Minimizing the function G(s ) we suppose that 

01+21 (/?2+"2)(P+ JV (P-3)) 

s = r(T) = G(s , m in) = c(0) maxlT^^ii+^ic- 1 ), T( i +/32)(P+ iv( P -3))+2Jv /32 } = 

1 p+iV(p-3) 

c(0) max{Tp+"(f>- 2 ), T^c- 2 )}. (3.6) 

Thus Uk(T,x) = for all x G {x = (x',X]sr) : ^ r(T)}, k — 1,N. And 
Theorem 1.2 is proved completely. 

□ 



Appendix A 

Lemma A.l. [T8] Le£ £/ie nonnegative continuous nonincreasing function 
f(s) : [so, oo) — > M 1 satisfies the following functional relation: 

f(s + f(s))^ef(s)Vs>s , 0<e<l. 

T/ien /(s) = Vs ^ s + (1 - e)" 1 /^). 

Lemma A. 2. [9\ If ft C is a bounded domain with piecewise- smooth 
boundary, a > 1, b G (0, o), and d > 1, then there exist positive constants 
d\ and d 2 (d 2 = if the domain Q is unbounded and d 2 = d 2 5 ^ if the 
domain Q = K(s, 5)) that depend only on Q, d, b, and N and are such that, 
for any function v(x) G W 1,d (Q) D L b (Q), the following inequality is true: 

IMIz/*(fi) ^ di \\Dv \\ L d(fy Iklli6(n) + d 2 \\v \\ L b^ 
where 9 = j b j_ a j G [0, 1). 

6 + lV~d 
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